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Abstract 



In this paper we investigate the hadronic mass spectra of inclusive B de- 
cays. Specifically, we study how an upper cut on the invariant mass spectrum, 
which is necessary to extract Vub, results in the breakdown of the standard 
perturbative expansion due to the existence of large infrared logs. We first 
show how the decay rate factorizes at the level of the double differential distri- 
bution. Then, we present closed form expressions for the resummed cut rate 
for the inclusive decays S — > X<j 7 and B — > Xu eu a.i next-to-leading order 
in the infrared logs. Using these results, we determine the range of cuts for 
which resummation is necessary, as well as the range for which the resummed 
expansion itself breaks down. We also use our results to extract the leading 
and next to leading infrared log contribution to the two loop differential rate. 
We find that for the phenomenologically interesting cut values, there is only 
a small region where the calculation is under control. Furthermore, the size 
of this region is sensitive to the parameter A. We discuss the viability of 
extracting Vub from the hadronic mass spectrum. 
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I. INTRODUCTION 



Inclusive B decays are considered fertile ground for precision tests of the standard model. 
The process B e v can be used to extract the all important Cabibbo-Kobayashi- 

Maskawa (CKM) matrix element Vub-, while B decays are important for discover- 

ing new physics. However, the utility of experimental measurements of these processes is 
bounded by our ability to control the theoretical errors. Tremendous effort has gone into 
determining ways to calculate these rates in a systematic fashion. Indeed, the algorithm 
for calculating these rates is now part of the theoretical canon Unfortunately, exper- 
imental cuts complicate life for theorists. In particular, these cuts often force us to work 
near the boundary of the phase space, where the aforementioned canonical techniques break 
down. Only now are we learning how to retool our calculations to accommodate these highly 
non-trivial issues. 

The complications arise if the cut forces us into a corner of phase space, since the calcula- 
tion can now depend on a new parameter, p, which is a measure of the relative size of phase 
space of interest. When this parameter becomes parametrically small, the systematics of the 
calculation usually break down; perturbative QCD corrections become enhanced by large 
logs of the form logp, while the non-perturbative expansion in A/m^ becomes an expansion 
in A/(pmb). 

The most relevant cut rate arises in semi-leptonic B decays, where one wishes to mea- 
sure Vub by eliminating the large background from charmed transitions. To eliminate this 
background, we have a choice of variables with which to cut. Perhaps the simplest choice is 
the electron energy, which is the oldest method used for extracting \Vub\- Unfortunately, as 
has been widely discussed in the literature, such extractions are typically model dependent 
since the rate in this window is sensitive to the Fermi motion of the heavy quark. There 
is no way to write down a meaningful theoretical error for such extractions. It is only 
very recently that a model independent method has been proposed, within a well defined 
systematic scheme 0, that could lead to an extraction with a well defined error. 



2 



It is also possible to remove the background from charmed transitions by cutting on the 
hadronic invariant mass While this choice presents a greater experimental challenge, it 
benefits from the fact that, unlike the electron spectrum, most of the B —>■ X^ev decays 
are expected to lie within the region sh < M"^. Furthermore, it is believed that even 
though both the invariant mass region sh < Mjj and electron energy regions Mb /2 > Ef,> 
(M^ — Mjj) / {2Mb) receive contributions from hadronic final states with invariant mass up 
to the cut mass spectrum will be less sensitive to local duality violations. This belief 
rests on the fact that the contribution of large mass states is kinematically suppressed for 
the electron energy spectrum in the region of interest. 

The goal of this paper is to study the viability of extracting Vub from the invariant mass 
spectrum. In particular, we are interested in studying the breakdown of the perturbative 
expansion for the cut rate, and whether or not a reorganized expansion can be used reliably. 

Building upon the work of Korchemsky and Sterman 0] , we begin by discussing how the 
doubly differential decay rate factorizes in moment space. We then use the recent results of 
to calculate a closed form expression for the inverse Mellin transform at next-to-leading 
logarithmic (NLL) order. The result for the resummed rate is presented in terms of the 
partonic as well as hadronic invariant mass. This result is used to extract a piece of the 
two loop rate, the size of which can be compared to the Brodsky-Lepage-Mackenzie (BLM) 
two loop correction. We then determine the region of invariant mass where resummation is 
necessary as well as the region where the reorganized expansion breaks down. We conclude 
with a brief discussion of the phenomenology, saving a complete discussion, including the 
effects of Fermi motion, for a later publication |^. 

II. FACTORIZATION IN HADRONIC VARIABLES 

The problem of summing large threshold logarithms in perturbative expansions, which 
arise due to incomplete KLN cancellation of the IR sensitivity at the edge of phase space, has 
been addressed for various processes [Q. The technique rehes on factorization, which allows 
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for resummations via a renormalization group equation. The factorization in B decays has 
been previously discussed in leptonic variables in [^,0. Here we review the arguments that 
are germane to our discussion of factorization in terms of hadronic variables. 

Consider the inclusive semi-leptonic decay of the b quark into a lepton pair with momenta 
Q = {Pe + Pu) and a hadronic jet of momenta ph- It is convenient to define the following 
partonic kinematic variables in the rest frame of the B meson v = (1, 0), 

Ph , "^v-ph 2vpe 

so = —, h= , x= , (1) 

rrib mi, 



with phase space boundaries 



h 1 r- — h 1 



0<so<l, 2^J rso < h < 1 + rso, l----V/i'-4so<x<l-- + -V/i2-4So. (2) 

In addition, it is customary to define the leptonic variables 

2v-q , , 

vo = , y = —- (3) 

TJib mi 

In terms of the leptonic variables, sq = {1 — + y) and h = 2 — y^, one can see that in the 
endpoint region of the electron energy spectrum when x —>■ 1 with y < 1, the invariant mass 
of the jet approaches zero with its energy held fixed. In addition, the jet hadronizes at a much 
later time in the rest frame of the B meson, due to the time dilation. Factorization exploits 
this and separates the particular differential rate under consideration into subprocesses with 
disparate scales. This factorization fails when the jet energy vanishes in the dangerous region 
y ^ X —>■ 1. However, this problematic region of phase space is suppressed because the rate 
to produce soft massless fermions vanishes at tree level. 

The infrared sensitive regions, which give rise to the large logarithms, can be determined 
by constructing a reduced diagram, as shown in Fig. 1. According the Coleman- Norton 
theorem ||^, a diagram at the infrared singular point must describe a physically realizable 
process after contracting all off-shell lines to a point. In the figure, 5* denotes a soft blob 
which interacts with the jet and the b quark via soft lines. J denotes the hadronic jet 
and H the hard scattering amplitude. Thus, the reduced diagram in Fig. 1 is simply a 
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FIG. 1. Reduced diagram for inclusive B decays. 



visualization of factorization. An important consequence of this factorization is that the soft 
function S is universal. Moreover, it has been shown [^] that the soft function S contains 
the non-perturbative structure function introduced in Ref. [^. It is this universality that 
will eventually allow us to eliminate the dependence on unknown non-perturbative hadronic 
dynamics @,||. 

From the above discussions one can see that in terms of the variables introduced earlier, 
factorization holds when sq ~ and /i ~ 1, or equivalently 1 ^ s^/h. The typical momenta 
flowing through the hard subprocesses are 0{mf,). Thus, H does not contain any large 
threshold logarithms and has a well-defined perturbative expansion in as{mh). The soft 
function S contains typical momentum k, with k~ ^ k± = 0{mhSQ/h). By soft we 

mean soft compared to mj,, but still larger than Aqcd- For an energetic u quark moving 
in the '— ' direction, the jet subprocess has typical momenta p such that p" ^ P^,P± with 
p^ = 0{mb sq/H), p]_ = 0{ml so/h), and p^ = 0{mh). In order to delineate between 
momentum regimes, a factorization scale /i is introduced. The fact that the process is 
independent of the factorization scale fi is utilized to sum the large threshold logarithms in 
the soft and jet functions. The reduced diagram for the inclusive radiative decays B —>■ X^j 
is exactly the same as above if we replace the lepton pair with a photon and ignore the 
strange quark mass. 

In terms of Sq and h, the triply differential rate, which factorizes into hard, jet and soft 
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subprocesses 0], may be written as 
1 d^T 



Tq dso dh dx 

fMB/mt 

12{2-h-x){x + h-l-So) dzS{z)m;j[mf,h{z-^),fi]H{mbh/fi), (4) 

ro = Y^IK.«, (5) 

where ^ = 1 — sq/Zi is analogous to the Bjorken scahng variable in deep inelastic scattering. 
z = l + k^/rrih, where is the heavy quark light cone residual momentum. S{z) essentially 
describes the probability for the b quark to carry light cone momentum fraction z and allows 
for a leakage past the partonic endpoint, as can be seen explicitly in the upper limit of z. A 
similar factorized expression holds for the inclusive radiative decays near the endpoint, 

1 dr^ fMB/m, 



To dso 



rlviB/'irib 

/ dzS{z)mlJ[ml{z-l + So),fj.^]H{mb/fJ.), (6) 

J l—sn 



G 



2 



2/1^(3 -2/i) / dzS{z)mlJ[mlh{z-0,fJ']H{mbh/fi) + 0{so). (8) 

J h 



^o = ^,\VM"^C^^l (7) 

For the inclusive semi-leptonic decays, the integration over x can be done in the endpoint 
region and the resulting doubly differential rate i^ 

Tq dso dh 

This is where the factorization in the invariant mass spectrum is simpler than in the electron 
energy spectrum. In the latter case, none of the integrals in the triply differential rate can 
be done trivially and one has to take an extra derivative with respect to x to arrive at an 
expression similar to Eq. An interesting consequence of factorization is that it connects 
the electron energy spectrum in the region x 1 with the invariant mass spectrum, in the 



region sq — > 0, 



dso 



1 d dT , . 

, 2 dx dx^ 



^ Here we confirm explicitly that the dangerous region /i — > 0, where the energy of the hadronic 
jet vanishes and factorization fails, is suppressed by the pre-factor and therefore not important. 
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which can be verified exphcitly at the one loop level using the corresponding expressions in 
Ref. dig. 



III. THE PERTURBATIVE RESUMMATION 



At one loop level the differential rates are |10|JTl 
1 d^T 



Tq dso dh 
1 dV^ 



2h' (3 - 2h) 



Fq' dso 



1 



Ait 



4 . 
13 + -vr^ 



+ — E2[h, So) 
47r 



CpCXs 

in 



7 



6 + 3so-2s^-2(2-so)logSo- — + 4 



.So 



l og -So ' 

^0 , 



(10) 



where 



Ei{h,so) = -5(so) 



8 \og\h) - 10 log h + + 4 Li2(l - /i) + 5 + ^vr^ 



81og(i±^)+7(l-t) 



rn /, N c-/-. N 4/l^lOg/l 

E2(/i,So) = 5(so)— ^-4 



1 - h 

+ 4/it(10- 15/i + 8so), 



2/i(3 - Ah) - 3(1 - 2h)sQ - 2s^ log 



1 +t 
l-t 



(12) 



(13) 



and t = Jl — Asq/K^. We also adopt the following definition for the '+' distributions 



|^ log"(go) ^ 

This definition is such that 
'log"(So)^ 



lim 



, log"(go) 
■So 



+ 5 (so; 



log' 



n+li 



n + l 



(14) 



fp 

/ dso-F(so) 

JO 



log"+Vp) /-p 
m -^-^ + / c^so [F{so) - F{0)] 
n + l Jo 



log" (so) 



(15) 



So / ^ n + i Jo Sq 

Note that if the parameter p becomes parametrically small, the first term on the right 
hand side of Eq. ( p!5[ ) will give large logarithms thereby spoiling the systematics of the 
perturbative expansions, while the second term must be regular as p — > 0. To perform the 
resummation we go into the moment space where the amplitudes factorize completely. In 
the case of inclusive semi-leptonic decays, it is convenient to define a new variable X = so/h 
with kinematic range 



Mi 



N 



Jo in 



The region with h > 1 is populated with real gluon emissions only, whereas the region with 
h < 1 has both real and virtual gluon corrections. In the relevant region h ^ 1 and A ^ 
the contributions from real and virtual gluon emissions combine to give terms which are '+' 
distributions, which upon integrating up to a cut, lead to the large logs we wish to resum. 

To proceed, we take the iVth moment with respect to ^ = 1 — A in the large N limit. In the 
region sq ~ and z ~ ,^ ~ 1, one can replace J[ml h{z — ^)] in Eq. (|D with J[ml h{l — C,/z)]. 
This replacement is permissible to the order we are working. We then obtain 

/o ^ ' Tq dX dh 
= 2/1^(3 - 2h) Sn JN{ml h/jj'^) H{mi, /fi) + 0{1/N), (17) 

Mml/fi^) = ml f'dyy''-'j[ml{l-y),fi], (18) 
Jo 

Sn= / dzz^Siz). (19) 
Jo 

The soft moment Sn- further decomposes into a perturbative soft piece, which accounts for 
soft gluon radiation and a non-perturbative piece which incorporates bound state dynamics 
and serves as the boundary condition for the renormalization group equation 0]. 

S{z)= -y-f[m,{l-y)]a{z/y), (20) 

Jz y 

Sn = Jn (7n, (21) 

where f{y) = {B{v) \ S{y—iD^) \B{v)) is the non-perturbative structure function defined 
in Ref. ||^. A similar expression holds for the inclusive radiative decays B Xgj 

/■I 1 dV^ 

Ml= drs,ii-sof-'^^ 

Jo i dso 

= JnctnJnH^ + 0{1/N), (22) 

Subsequently, we will ignore the non-perturbative structure function f{y) and concentrate 
on the perturbative resummations. 
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It merits emphasizing that the large asymptotics of the moments corresponds to the 
behavior of the spectra in the region sq ~ A ~ 0. Taking the large limit also enables us 
to extend the integration limit of A in Eq. (|1^) up to 1, despite the fact that the kinematic 
range of A never goes up to 1. In this limit the contribution from the region A ~ 1 is power 
suppressed. 

Comparing Eq. ([T7|) with the corresponding expression for the electron energy spectrum 
in Ref. one sees that the moments (Tn and Jn are identical with those in the electron 
energy spectrum, with change of variables x — > 1 — A and 2 — yo^h. A similar identification 
for the resummed radiative decays, Eq. ([2^), can be made with the change of variables 



X ^ 1 — Sq. In moment space the soft and jet functions [|T2| have been calculated to NLL 
order and are given by 

aN Jn = exp [ \og{N) g, (x) + 92 ix)] (23) 

Jl = exp[log(iV) g,ix) + ^^^(x)], (24) 

where x = «s('^b)/5ologiV, and gi and g2 are given explicitly as 0] 

9iix) = -TT^Kl - 2X) log(l - 2x) - 2(1 - x) log(l - x)], (25) 

92{x) = glix) + 9si{x,h), (26) 



1 



log(l-2x)-21og(l-x) 



+ ^log2(l-2x)-log2(l-x) 



1 2 

iog(i - x) - iriT iog(i - 2x) 



+ pf- [log(l - 2x) - log(l - x)] , (27) 
9 Ax, h) = :r^ log(/i) log(l - x). (28) 

In the above, (3o = (11 Ca - 2 iV/)/(127r), /3i = (17 - SC^iV/ - 3Cf Nf)/{2ATT^), k = 
Ca (67/18 - 7rV6) - lOT^j A^//9, and -/e = 0.577216... is the Euler-Mascheroni constant. In 
our case, Ca = 3, Cp = 4/3, and Tr = 1/2. To get back the physical spectra from the 
moment space, the inverse Mellin transform has to be evaluated at NLL accuracy as well. 
To this end, we apply the identity derived in the Appendix of Ref. [@ 
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1 /'C+200 

— / ^-N ^\og{N)Fi[as\ogN]+F2[as\ogN] 

2TTi Jc-ioo 

where / = — log(— logx) ~ — log(l — x\ and 

oo fe— 1 

^K/) = 5:«^^^p (30) 

k=\ j=0 

represents next-to-next-to-leading log contributions. Changing variables from A back to 
So/h, we obtain 
1 d^T 



2h\?>-2h)H{h) 



Tq dso dh 

dso\ \h r[l - gi{a,f3ol) -asPolg[{as(3ol)]f ' 

1 dY'r d f „ gi5i(a./3on)+sJ(a./3oP) ^ 

f^^^-^^dfoX ^^'^ ~ r [1 - g,{as(3o I'') - «s/?o g'MsPo l^)] J ' ^^'^^ 
where / = — log (— log(l — A)) ~ — log(so//i), and V = — log (— log(l — sq)) ~ — log(so). 
The ^-functions define the differential rates in a distribution sense, as ?7 ^ 0, and turn the 
singular terms into the '+' distributions, as can be seen explicitly by expanding in power 
series of log sq and using the definition Eq. (|T^) . 

The hard parts can be obtained through the one loop results Eq. (|10D and Eq. ( pT]) 



Uog\h) - 51og/i + + 2Li2(l -h) + - + — 

B K J ^ 3-2h ^ '23 



(33) 



Eq. ( pTD and Eq. ( |32D reproduce the dominate contribution at one loop level in the limit 
A ^ and sq — > 0, respectively, and include the infinite set of terms of the form a" log""*"^(so) 
and a;"log"'(so) in the Sudakov exponent for both semi-leptonic and radiative decays. 



IV. THE INTEGRATED CUT INVARIANT MASS SPECTRUM 



As previously mentioned, it has been proposed that we measure the modulus of the 
CKM matrix element Vub from inclusive decays by making a cut on the hadronic invariant 
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mass below M|,. This cut eliminates the overwhelming background from bottom to charm 
transitions. While it is the hadronic invariant mass which is of interest, we shall first consider 
the cut partonic invariant mass, as it will be relevant to our conclusions. 

The use of an upper cut cq on the partonic invariant mass introduces large logs of the 
form as log^(co). As Cq approaches zero, the logs become parametrically large and need to be 
resummed. Using the resummation formulas in the previous section, it is simple to generate 
an expression for the resummed cut rate, since our expression can be written as a total 
derivative with respect to sq. The cut rate may be written as 

I- r(co) = £ dh 2h\3 - 2h) [e (I) - U (I)] + dso dh 7„.(/i, so), (35) 



r 

where 



gi(as/3o 0+ff2(as/3o 



(36) 

l=- log(co//i) 



r [1 - giiosPo I) - asPo I g'l {asPo I)] 
and 7a^(/?., So) is simply the one loop rate defined in Eq. ([T0|). stands for terms up to 

0{as) when expanding in power series of as- We subtracted it from ^ in order to ensure 
that we correctly reproduce the one loop result at order a^. From Eq. ( pS] ) one can see that 
the large logs arise when Sq approaches zero, which is not only the lower kinematic limit for 
So, but also the phase space boundary for virtual and real gluon emissions. 

The perturbative expansion has been reorganized into an expansion in the exponent. 
The systematics of this expansion have been discussed at length in PJT3[|. Here we just 
recall that we may test the convergence of the reorganized expansion by comparing the NLL 
resummed result with the leading-logarithm (LL) resummed result. Fig. 2 shows the cut 
rate as a function of cq. In this figure we show the one loop cut rate as well as the resummed 
cut rate with and without NLL corrections. We see that for 0.1 < Co < 0.2 the resummation 
becomes necessary, while for cq < 0.1 the NLL dominates the LL, so that we can no longer 
trust our results. Indeed, this breakdown occurs well before we reach the Landau pole at 
Sp = e-V(2a./3o) ~ 0.028. 

We now expand this result to pick out the leading and next to leading infrared log 
contribution to the two loop differential rate. This contribution is given by 
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FIG. 2. The rate as a function of the partonic cut. The dotted line is the one loop result, the 
dashed line is the LL result, while the solid line is the NLL log result. The difference between the one 
loop and resummed results at large cq is due to two loop corrections introduced in the resummation. 
We use as{rnb) = 0.21. 



dr 



Tq dsodh 



2h'^{3-2h) 



37r2 



8 /log=^(so) 



+ (14 - GnPo) 



(37) 



As expected, we see that the most singular contribution at 0(«^) doesn't have any terms 
proportional to (3q. It may be the case that, in this particular region of phase space the 
infrared logs terms may dominate over the BLM terms. Such a conclusion was reached 
in [p|,[l3| for the two loop contribution to lepton and photon spectra in semi-leptonic and 
radiative decays, respectively. Naturally , this does not preclude the possibility that there 
exist a cancellation with other uncalculated terms, such that the a^jSo still dominate. 

Let us now consider the physical case, where we are interested in placing an upper cut 
on the hadronic invariant mass. The hadronic invariant mass may be written as 



sh 



sh 



so + eh + 



(38) 
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where A is the mass difference Mb — mi, in the infinite b quark mass hmit, which is a measure 
of binding energy for the b quark inside the B meson. Thus, given a cut on sh, c, we may 
translate this into an h dependent cut on sq. After changing the order of integration we find 
that the cut rate may be written as 




if c > e(l + e), whereas if c < e(l + e) then 




This situation is different from the partonic invariant mass spectrum. In the present 
case, the lower kinematic limit of sh is e^, while the phase space boundary for virtual 
and real gluon emissions has been moved to = e(l + e), above which only real gluon 
emissions contribute.^ Hence there could potentially be two kinds of parametrically large 
logarithms log(c — e^) and log(c — pe) appearing. However, as can be seen from Eq. (p8|), 
Sh corresponds to vanishing sq and h, which is the dangerous region where the infra- 

red factorization fails and is kinematically suppressed by the tree level rate. More explicitly, 
log(c — e^) will be killed by the pre-factor 2(c — e^)^ so that the rate vanishes as c — e^. 
The only important parametrically large logarithms are those resulting from the incomplete 
cancellation between virtual and real corrections at the phase space boundary p^. From this 
argument it would seem that we only need to resum logs of the form log(c — p^). However, in 
this particular case, the distance from p^ to is only e, which is expected to be a numerically 
small quantity at around 0.08. If the experimental cut c lies below p^, the partially integrated 
rate could still be sensitive to log(c — e^) due to the smallness of e. We thus chose to resum 



more general discussion for this kind of phenomenon in other observables can be found in 
Ref. [|14|. 
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all logs of the form log(c — e^) as well. Fortunately, since we have the resummed rate at the 
doubly differential level, this is not a problem. The resummed rate with hadronic mass cut 
c> pe is given by 

1 



-Tic) 



where 







dh2h'{3-2h) U 



c - 



h{l + e) 



+ / ds 



e — Pe 

1 + E 



/ dh2h\3-2h) ^'{so/h)-C^X^o/h) 



c—pc 
1 + e 



+ / dsQ f dh'yaXh,So)+ I ds 
/o 



1 + e 



/ d/i7„^(/i. So), (42) 

2v^ 



(43) 



While the resummed rate with hadronic mass cut c < is given by 



^r(c) 

J- 



dh2h^{3-2h) U 



h 



+ dso lasihJo) 

Jo J^yT^i 



(44) 



Analytic expressions for the partially integrated rate at one loop level, last lines in Eq. (^ 
and Eq. (H), can be found in Ref. [TO . 



V. RESULTS AND CONCLUSIONS 

In Fig. 3 we show the one loop, LL resummed (including gi only) and NLL resummed 
(including gi and 5^2) results for A = 0.39 GeV. We see the for c < 0.18, the next to leading 
order approximation breaks down, as the next to leading order piece becomes just as large 
as the leading order piece. Notice that when c ~ 0.18, the effective cut on the partonic 
invariant mass is cq — 0.09, which from Fig. 2 we see is consistent with the breakdown of 
the resummed expansion. This is contrast with the result that an energy cut of 2.1 GeV, 
on the rate for B Xs'y does not necessitate resummation, as the argument of the logs in 
the case of the radiative decay is about 0.12. 
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FIG. 3. The rate with a hadronic cut. The dotted line is the one loop result, the dashed line is 
the LL resummed result, while the solid line is the NLL resummed result. For A = 0.39 GeV, we 
have pe = 0.0861. In the region above p^, we run into Landau pole at Cp = 0.1159 and interpolate 
in the small region between /9e and Cp. 

It is clear from Fig. 3 that resummation shifts the whole spectrum toward the high 
invariant mass region such that the number of events which lie below the cut c is decreased. 
This occurs because the high invariant mass region with sh > Pe is populated with real 
gluon emissions only. 

In Fig. 4 we show the cut rate for several different values of A. It is clear that the cut 
rate can be very sensitive to the value of the unphysical parameter A. This occurs because 
the argument of the large logs is now c — p^. The parameter A is well defined at a given 
fixed order in perturbation theory and the values chosen in Fig. 3 are the mean value and 



one sigma values extracted in Ref. [15] at one loop. Soon the errors in the extraction of 
A will become less significant. In this paper we will not delve into the phenomenology 
of the extraction of V^b, as the calculation we have discussed here has not included the 
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FIG. 4. The cut rate for several different values of A. The dashed line is for A = 0.28 GeV, the 
dot-dashed line is for A = 0.50 GeV, while the solid line is for A = 0.39 GeV. The dotted line is 
the one loop result for A = 0.39 GeV. 

important non-perturbative corrections coming from the Fermi motion. These corrections 
are parameterized in terms of a well defined structure function, f{y) defined in Eq. (pOD, 
and as discussed in [0], can be very large. We will reexamine the issue of the structure 
function in a future publication, where following , we will eliminate the structure function 
from the cut rate prediction by utilizing the data from the end-point of radiative B decays, 
which in turn encodes all the information contained in the structure function. 

From Fig. 4 we see that, for large values of A, the resummation is not under control for 
phenomenologically interesting cuts and thus it may not be possible to extract Vub in this 
way. However, when modding out by the soft function using the b ^ s'y rate, there are 
cancellations which lead to a perturbative series which is better behaved. This is indeed 
what happened in the case for the electron energy spectrum 0. Thus, we refrain from 
drawing any conclusions regarding the viability of extracting Vub from the invariant mass 
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spectrum at this time. 

The purpose of this paper was to determine the effect of threshold resummation on the 
rate for semi-leptonic B decays with a cut on the hadronic invariant mass. We first showed 
how the rate factorizes when written in terms of hadronic variables, generalizing the results 
of 1^. Using this factorization, we resummed the cut rate at next to leading order in the 
infrared logs and found that, for cuts of interest, the resummation is crucial, and that for 
c < 0.18, even the next to leading order resummed rate is no longer reliable. However, this 
breakdown point depends on the value of A, and becomes smaller as A is decreased. A more 
phenomenological analysis, including the effects of the structure function responsible for the 
Fermi motion, is forthcoming. 
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